Surface Waves on a Semi-toroidal Water Ring 



Sunghwan Jung ^, Erica Kim ^, Michael J. Shelley ^, and Jun Zhang ^'^ 
^ Applied Mathematics Laboratory, Courant Institute of Mathematical Sciences, 
New York University, 251 Mercer Street, New York, New York 10012, USA 
^ Department of Physics, New York University, 
4 Washington Place, New York, New York 10003, USA 

We study the dynamics of surface waves on a semi-toroidal ring of water that is 
excited by vertical vibration. We create this specific fluid volume by patterning a 
glass plate with a hydrophobic coating, which confines the fiuid to a precise geo- 
metric region. To excite the system, the supporting plate is vibrated up and down, 
thus accelerating and decelerating the fluid ring along its toroidal axis. When the 
driving acceleration is sufficiently high, the surface develops a standing wave, and 
at yet larger accelerations, a traveling wave emerges. We also explore frequency 
dependencies and other geometric shapes of confinement. 

In 1831, Faraday first observed that surface waves on a vibrated fluid volume oscillate 
at half the driving frequency [l|. These Faraday waves demonstrate the surface instability 
caused by strong oscillatory accelerations exerted on the fluid volume. Benjamin & Ursell 
showed that the amplitude of a surface eigenmode obeys the Mathieu equation Q3. It 
follows that the superposition of multiples of harmonics and subharmonics are also solutions 
to the Mathieu equation. In the presence of finite viscosity, however, the subharmonic 



response is dominant 



Previous experiments have revealed that patterns of various symmetries, such as striped 
0,3, triangular S, square Q, fl \\\ , and hexagonal [is], [l^ , can be excited on the free 
surface of a fluid layer. Quasi-one dimensional surface waves have also been studied in both 
narrow annular and channel geometries Q, 17, 18]. In these experiments, the standing 
waves often interact with the meniscus that forms at the bounding walls. The contact point 
and length-scale of this meniscus continually change as a result of the vertical oscillations. 
Consequently, the meniscus emits waves towards the bulk [l^. In other experiments, to 



remove this meniscus effect, the contact point is pinned on a sharp edge or brim 
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2J, |25|, l26|, has been 



A curved fluid surface, such as on a hemi-spherical drop l22l, l23 
used to study surface waves in a conflned geometry without bounding walls. Vibrating the 
drop causes waves to form on the curved surface, and for high external forcing, droplets are 
ejected. This process, termed atomization, provides one way to create a spray 27|, l28|]. The 
instability of the drop's semi-spherical surface has been utilized for spray cooling, mixing, 
and humidiflcation. 
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27l |. The differences be- 



Faraday instability has been well studied both in flat fluid sur 
by the lateral boundaries and in sessile (pendant) drops |23|, l24j . 
tween the two cases are their aspect ratios and boundary conditions (moving lateral contact 
lines or the capillary pinning). The geometry of our interest lies between the two cases: the 
semi-toroidal geometry is spatially extended with high aspect ratios (system size/capillary 
length) and pinning boundary conditions apply. It is however much less studied compared to 
the above two settings. Here, we employ a hydrophobic/hydrophilic patterning of a surface 
to conflne a water volume to a speciflc region. In particular, if part of a hydrophilic sub- 
strate (low surface energy) is covered by a hydrophobic coating (high surface energy; we use 
Fluorothane ME; Cytonix Corp.), water is prevented from spreading beyond the hydrophilic 
region. The boundary between the two regions acts as an edge upon which the fluid contact 
point is pinned. 

In the experiment we report here, the hydrophilic region is a thin annulus, and water 
placed there forms into a semi-toroidal volume with two pinned contact lines. We flnd 
that surface waves flrst develop along the centerline of the torus, which allows us to study 
quasi-one dimensional waves on a curved surface. 

The schematic of our experimental setup is shown in Fig. [H The annular hydrophilic 
region has inner radius ai = 2.5 cm and outer radius 02 = 3.5 cm. The patterned glass plate 
is rigidly connected to a speaker, which oscillates the plate vertically at a controlled driving 
frequency and amplitude. A high-speed video camera is centered above the plate along the 
axis of the annulus. Since the semi-torus of water has a cross-sectional radius (r ^ 0.5 cm) 
that is relatively small in comparison to the distance between the plate and the camera, 
bright regions of the water ring seen in top-viewed video frames correspond to local extrema 
of the fluid surface height. 

We use a systematic procedure to deflne the onset of surface waves. At flxed driving 
frequencies Q/27r between 20 Hz and 65 Hz, the amplitude of oscillation is slowly increased 
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FIG. 1: (a) Side-view of the experimental set-up: a semi-toroidal water-ring on the hydrophobic 
patterned glass plate that is vibrated vertically by a loudspeaker, (b) Top and cross-sectional views 
of the water-ring. 

until the fluid surface loses stability to azimuthally modulated standing waves. Distinct 
wave patterns are observed at half the driving frequencies. Figures [2^ and b show these 
standing waves for two different driving frequencies. At 22 Hz (a), the wave pattern has 12 
nodes (that is, 12 crossings of the wave peak along the annular mid-line). The left half of 
Fig.[2](a) shows a snapshot of the surface wave, whereas the right half is the image of two 
superimposed surface waves in different phases. In the snapshot, the cross- sectional shapes, 
or onset modes, of the fluid surface are shown as insets when cutting across the I and II 
lines. The overlapped image demonstrates that standing waves are formed at half the driving 
frequency (22 Hz) near the onset. Increasing the driving frequency to 36 Hz increases the 
number of wave nodes to 24 (see Fig.[2]D). Near the onset, as the driving frequency increases, 
the number of nodes increases monotonically (and discretely, due to the geometry), as is 
shown in Fig. [3l 

To describe the instability to the azimuthally modulated surface waves seen in Fig. [21 we 
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FIG. 2: Images of surface waves emerging from the water ring. At the driving frequency 22 Hz, 
12 nodes are excited in the azimuthal direction (a), whereas 24 nodes appear at 36 Hz (b). In 
(a), the cross-sectional shapes along lines I (referred to as the peak) and II (trough) are shown 
schematically. The dashed curves represent the fluid surfaces in equilibrium. White dots represent 
the cross-section position of the bright curve in a top-viewed image. The right half of image (a) is 
an overlap of two photographs differing by the period of excitation T (=1/22 sec). 

consider an approximate model which captures the main aspects of our observations. First, 
we neglect viscous effects. Second, instead of using the exact semi-toroidal geometry for 
our calculations, we neglect the annular curvature and use a cylindrical geometry. Finally, 
we note that the capillary length (-^ 2 cr/p^') in our experiment is ~ 0.3 cm slightly smaller 
than the radius r ~ 0.5 cm, and so we expect the base state (i.e. the surface profile without 
azimuthal modulation) to have a flattened, temporally oscillating surface profile. Solving for 
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FIG. 3: Number of nodes around the annulus observed near the onset of wave formation. Dots are 
from experiments. The solid line is a result of the inviscid coupled linear model and the dashed 
line is from the classical dispersion relation (($7/2)^ = gk + ak^/p). 




FIG. 4: Lateral amplitude of oscillation vs. non-dimensional acceleration {hVL^ /g) at r2/27r = 30 
Hz. A is the transverse amplitude as indicated in the inset. This response curve shows hysteresis 
over driving acceleration between 0.63 and 0.67 g. 

the temporal evolution of the surface profile, even azimuthally unmodulated, is a nonlinear 
problem. Hence, for simplicity, we assume that the base state is close to a semi-circle with 
constant radius R. Small fluctuations of the surface can then be expressed as R-\-C,{9), where 
C,{9) is a small radial displacement deviating from R. This assumption leads to a simple 
linear analysis for the perturbed fluid surface. A fully nonlinear, numerical treatment of the 
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oscillations and Faraday instabilities of a driven viscous drop has been considered by James 
et al ui- 

A linearized form of Bernoulli equation with a velocity potential p] in a moving frame 



IS 



gXcosO +^--^^ + ^^ + ^J=0, (1) 
VV = . (2) 

Here, ge{t) = g + acosQt is the oscillating body force where g and a are the acceleration 
due to a gravity and the peak acceleration of the vibrating plate, respectively. By using 
Vr\r=R = dr4>\r=R = dtC which ueglects the nonlinear terms in the kinematic boundary 
condition, we can assume solutions such as 



C{e,z:t) = ^Uili)^**''' (3) 
^(..,.;«)^^^^e.<--.. (4) 

where k is the axial wavenumber (= 27i/L; L is the wavelength) and m is the azimuthal 
wavenumber (in the direction of 9). Here, we assume that a/g is small and neglect higher 
azimuthal modes (|m| > 2) since we observed that m = 2 is the dominate mode near the 
onset of surface excitation. Using the symmetry (m,k = —C-m,k, setting the unphysical 
volume change term (^ofc to be zero, and plugging Eqs. (Pfl|) into Eq. ([T]) yields 



Zikdtt + cTik ge{t) 

ge{t) Z2kdtt + CT2k 



)l 


( Cik\ 




\ C2k / 



0, (5) 



where Z^k = Im{kR)/kI'^{kR) and amk = (cr/p) ■ + rri^ / - l/R^). The surface 
deformation is assumed to be of the Floquet form {C,mk = e"^* ^ C^'^^e*"'^*) where r = s + iaQ. 



fl. 



For subharmonic solutions {a = 1/2), one gets the recursion relation: 

HnCt^ = ^ (C2 + CS!) + O {{a/gf) , (6) 
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where 

Hn = [-1 + 4(Zu,^^ _ a,k)iZ2kUjl - <y2k)/g^] (7) 

and uJn = (1/2 + n)fl. The problem of stability of the cylindrical fluid surface is reduced to 
finding the eigenvalues of Eq. ([6]) up to the order of a/g. 
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FIG. 5: Number of nodes vs. non-dimensional acceleration at n/2TT = 30 Hz. As the acceleration 
is slowly increased, the number of nodes changes from 20 to 22 when the acceleration is about 
1.71 g. When the acceleration decreases, the number of nodes returns back to 20 at 1.48 5. The 
states shown in triangles (top branch) are associated with the traveling waves. 

FigureO compares the observed patterns with the inviscid coupled linear model and the 
classical dispersion relation ( (^^/2)^ = gk + ak^/p ). Our coupled model is in better 
agreement with our experimental observation than the classical dispersion relation. 

We also find experimentally that the onset of surface waves is hysteretic. Figure H] shows 
the maximal displacement A in the horizontal direction between standing wave peaks over 
the course of an oscillation, as a function of dimensionless driving acceleration hQ'^/g, where 
h is the amplitude of plate oscillation. Crossing a critical amplitude from below shows 
the sharp transition to standing wave. This bifurcation is subcritical, and at intermediate 
accelerations {0.63 g < hD? < 0.67(7) two stable states are observed. 

Figure [5] shows yet another hysteretic transition in the number of surface nodes, at a 
higher dimensionless acceleration. Bistable surface waves with different numbers of nodes 
co-exist between 1.48(7 aiid 1.71 g. We find that the surface waves above this transition 
still oscillate at half the driving frequency, but now have more nodes and a finite drift 
speed (triangles). The inset images show the patterns of peaks (bright) and troughs (dark) 
in space-time. The arrows indicate the direction of azimuthal drift and its speed. If the 
acceleration is further increased to around 4(7, the traveling surface waves become more 
spatially localized, and droplets are ejected. 
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FIG. 6: (a) Below threshold, stable fluid surface constrained to three regions which are bordered by 
a hydrophobic coating, (b) Above the onset, unstable fluid surface vibrating at half of the driving 
frequency. At 50 Hz, fluid surfaces constrained to the shapes "A","M","L" become unstable at 
accelerations 2.51, 2.23 and 2.64 g, respectively. 

Our masking technique can be used effectively to make complicated hydrophilic regions, 
and to investigate the dynamics of fluid volumes trapped within them. Figure [6] shows 
three separate fluid volumes bound within domains (the letter shapes A, M and L) of dif- 
ferent shapes at both stable and unstable states. At fixed volume to base-area ratio of 0.16 
cm'^/cm^, which gives a fixed height for all shapes slightly above 0.18 cm, instability occurs 
at different thresholds for the different base shapes. The shape that with high aspect ratio 
(e.g. letter "L") appears to have the highest threshold to instability. The top image shows 
these volumes at rest, where the transparent glass surface of the hydrophilic region reveals 
a square grid beneath. The bottom image shows the surface wave dynamics of these fluid 
volumes as revealed, at least in part, by the distortions of the underlying grid lines. 

In this paper, we have studied the Faraday instability of a semi-toroidal water volume 
whose surface is initially flat in the azimuthal direction but tightly conflned in the radial 
direction. To constrain the fluid into a tight geometric area (an annulus), the base plate 
is patterned with hydrophobic and hydrophilic regions. This particular geometry allows 
us to study the surface waves arising from a quasi-one-dimensional system with a periodic 
boundary condition. Systems with high aspect ratios (system size/capillary length), as 
used in our experiment, have not been well studied previously. Despite of the high aspect 
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ratio, we find similar dynamics (such as the hysteretic excitation of surface waves) as found 
in systems with fiat surfaces. The classical theory is modified to describe the observed 
dispersion relation in our system. The result agrees with our experimental observations. 
For large forcings, we observe another transition, again hysteretic, to traveling waves with 
an increased number of nodes. Such a transition has yet to be understood theoretically and 
is a subject of future investigation. 
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